Interference effects in the photorecombination of argonlike Sc 3+ ions: 
Storage-ring experiment and theory 
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Absolute total electron- ion recombination rate coefficients of argonlike Sc 3+ (3s 2 3p 6 ) ions have 
been measured for relative energies between electrons and ions ranging from to 45 eV. This 
energy range comprises all dielectronic recombination resonances attached to 3p— >3d and 3p— >4s 
excitations. A broad resonance with an experimental width of 0.89 ± 0.07 eV due to the 3p 5 3d 2 2 F 
intermediate state is found at 12.31 ±0.03 eV with a small experimental evidence for an asymmetric 
line shape. From R-Matrix and perturbative calculations we infer that the asymmetric line shape 
may not only be due to quantum mechanical interference between direct and resonant recombination 
channels as predicted by Gorczyca et al. [Phys. Rev. A 56, 4742 (1997)], but may partly also be 
due to the interaction with an adjacent overlapping DR resonance of the same symmetry. The 
overall agreement between theory and experiment is poor. Differences between our experimental 
and our theoretical resonance positions are as large as 1.4 eV. This illustrates the difficulty to 
accurately describe the structure of an atomic system with an open 3d-shell with state-of-the-art 
theoretical methods. Furthermore, we find that a relativistic theoretical treatment of the system 
under study is mandatory since the existence of experimentally observed strong 3p 5 3d 2 2D and 
3p 5 3d4s 2 D resonances can only be explained when calculations beyond LS-coupling are carried out. 

PACS numbers: 34.80.Kw, 34.80Lx 



I. INTRODUCTION 

Dielectronic recombination (DR) 

e" + A q+ -» [A^-^ + ]** -► A {q -^+ + hv (1) 

is a two-step electron-ion collision process, where in a first 
step a multiply-excited intermediate state is created by a 
resonant dielectronic capture (DC, inverse Auger) and in 
a second step that intermediate state decays by photon 
emission. Another possible recombination mechanism is 
radiative recombination (RR) 



+ A q+ -> + hi 



(2) 



which by emission of a photon directly proceeds to a non- 
autoionizing state of the recombined ion. When both the 
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initial and the final state including the photons emit- 
ted are the same, RR and DR are indistinguishable and 
they may interfere, with the signature of interference be- 
ing an asymmetric DR resonance line profile in the en- 
ergy differential recombination cross section [jj], |^, ||, Q. 
Asymmetric line shapes have experimentally so far been 
observed only in partial recombination cross sections of 
highly charged uranium ions Such an observation 
is interesting from a fundamental point of view since it 
questions the widely used independent processes approx- 
imation (IPA) in the calculation of RR+DR recombina- 
tion cross sections. 



One motivation for investigating the recombination of 
argonlike Sc 3+ (3s 2 3p 6 ) ions with free electrons is the the- 
oretical predication by Gorczyca et al. j(| of a strongly 
asymmetric line shape for the 3p 5 3d 2 2 F DR resonance 
caused by interference between the recombination path- 



2 



ways 

e - +Sc 3+ (3s 2 3p 6 ) -» Sc 2+ (3s 2 3p 5 3d 2 ) 

\ I (3) 
Sc 2+ (3s 2 3p 6 3d) +hu. 

Theoretically they predicted a Sc 2+ (3s 2 3p 5 3d 2 2 F) reso- 
nance position of 14.6 eV and a width of 1.8 eV. The large 
width of the resonance is explained by the fact that it de- 
cays predominantly via a very fast Super-Coster-Kronig 
transition (reverse of the horizontal arrow in Eq. ||) with 
a calculated rate of 2.84 x 10 15 s^ 1 . 

A first attempt to experimentally determine the line 
shape of the Sc 2+ (3p 5 3d 2 2 F) resonance by measuring the 
total Sc 3+ recombination rate coefficient at a heavy-ion 
storage ring failed due to the rather large statistical un- 
certainty of that measurement |j (see also Sec. [il]). How- 
ever, distinct differences between theoretical and experi- 
mental peak positions were found and attributed to the 
fact that a correct theoretical description of the highly 
correlated 3p 5 3d 2 configurations is extremely difficult. 
Using Cowan's || atomic structure code we find that the 
Sc 2+ (3p 5 3d 2 ) configuration gives rise to 45 fine struc- 
ture components distributed over an energy interval of 
about 20 eV width. In order to account for correlation 
effects Hansen and Quinet Q considered the interaction 
of 10 initial and 16 final configurations in the calcula- 
tion of 3p 6 3d— ► 3p 5 3d 2 transition energies in Ca + . Still 
their result for the Ca+(3p 6 3d 2 D -> 3p 5 3d 2 2 F) transi- 
tion energy deviates 0.76 eV from the experimental result 
(29.34 eV) determined by photoionization of metastable 
Ca + (3p 6 3d) ions Jl§. This again illustrates the difficulty 
to accurately describe inner shell transitions in atomic or 
ionic systems with open 3d-subshells. In this situation 
DR measurements (and complementary photoionization 
measurements) can yield valuable spectroscopic informa- 
tion on such systems. 

Here we present new experimental as well as theoretical 
results for the total Sc 3+ recombination rate coefficient. 
The experiment is a twofold extension of our previous 
work MI. The new data have much smaller statistical 
uncertainties and they are extended to the much wider 
energy range 0-45 eV (previous range 12.2-18.2 eV) in- 
cluding all DR series limits due to 3p — > 3d and 3p — > 4s 
core excitations. 

We have also carried out both nonperturbative and 
perturbative calculations for the e~ + Sc 3+ photore- 
combination. The nonperturbative method is based on 
the rigorous continuum-bound transition theory and the 
close- coupling R- matrix approach ]Tl| . The perturbative 
evaluation is a treatment including radiative recombina- 
tion (RR), dielectronic recombination (DR) and their in- 
terference. Our calculations reveal that three adjacent 
resonances mask the low energy side of the 3p 5 3d 2 2 F 
resonance. This makes the experimental observation of 
the predicted |3| interference between the 3p 5 3d 2 ( 3 F) 2 F 
DR resonance and the continuous RR cross section dif- 
ficult. Moreover, we find that several experimental DR 



resonances not reproduced by the LS-coupling R-matrix 
calculations of Gorczyca et al. are due to relativistic 
effects. 

In Sec. [H] and Sec. HI we outline our experimental and 
theoretical methods. Experimental and theo retic al re- 
sults are presented a nd dis cussed in Sees. IV A and EVE , 
respectively. In Sec. IV C] possible origins for the asym- 
metry of the 3p 5 3d 2 ( 3 FpF DR resonance are explored 
by theoretically considering interference between DR and 
RR as well as interference between overlapping DR reso- 
nances. Our main results are summarized in Sec. [y|. 



II. EXPERIMENT 

The experiment was performed at the heavy ion stor- 
age ring TSR of the Max Planck Institut fur Kern- 
physik (MPIK) in Heidelberg. Details of the experimen- 
tal setup and the data reduction procedures can be found 

in Refs. @ || IB IJ . The 45Sc3+ i0nS W6re accelerated 
by the MPIK tandem booster facility to their final energy 

and subsequently injected into the TSR. The rather low 

charge to mass ratio of q/A — 1/15 leads to the rather 

unfavorable condition (see below) that the ion energy 

is limited by the maximum bending power of the TSR 

dipole magnets. With the maximum magnetic rigidity of 

Bp = 1.4 Tm the highest accessible ion energy can be 

calculated from 



E-JAk 931.5 MeV/u 




Bp 



A 3.107 Tm 



(4) 



to Ei/ A w 420 keV/u. In the storage ring the circulating 
18.9 MeV Sc 3+ ion beam was merged with the magnet- 
ically guided electron beam of the electron cooler. Dur- 
ing electron-cooling, the electrons have to move with the 
same average velocity as the ions. This condition defines 
the cooling energy E c = (E-J A)(m c /u) as 230 eV. In or- 
der to increase the phase space density of the already 
stored ions and thereby free phase space for the next in- 
jection pulse from the accelerator, electron-cooling was 
allowed to be effective for 2 seconds after each injection 
of ions. These 'ecool-stacking' cycles Q were repeated 
three times until the accumulated ion current stabilized 
at values of up to 5 /iA. Before data taking was started 
a prolonged cooling interval of 5 s after the last injec- 
tion allowed the ion beam to shrink to a final diameter 
of about 2 mm as verified with a beam profile monitor 
based on residual gas ionization [fj"7| . 

During the measurement the electron energy E e in the 
laboratory frame was stepped through a preset range of 
values different from E c thus introducing non-zero rela- 
tive energies 



E r , 



E c — v E c 



(5) 



between ions and electrons. In the data analysis a rel- 
ativistically correct expression p4l was used instead of 
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Eq. H Recombined Sc 2+ ions were counted with 100 +°% 
efficiency as a function of cooler voltage on a single par- 
ticle detector Q located behind the first dipole magnet 
downstream of the electron cooler. The dipole magnet 
bends the circulating Sc 3+ ion beam onto a closed orbit 
and separates the recombined Sc 2+ ions from that orbit. 
In between two measurement steps the cooler voltage was 
first set back to the cooling value in order to maintain 
the ion beam quality and then set to a reference value 
at E rc \ — 45 eV which is chosen to lead to a relative 
velocity where the electron-ion recombination coefficient 
is only due to a negligible RR rate. Under this condi- 
tion the recombination rate measured at the reference 
point monitors the background signal. Choosing short 
time intervals of 10 ms duration for dwelling on the mea- 
surement, cooling and reference voltages ensured that the 
experimental environment did not change significantly in 
between signal and background measurement. An addi- 
tional interval of 1.5 ms after each change of the cooler 
voltage allowed the power supplies to reach the preset 
values before data taking was started. 



At the rather low ion energy of 420 keV/u the cross 
section for electron capture from residual gas molecules 
is estimated — taking the measured residual gas com- 
position into account and using a semi-empirical formula 
p9[ for the charge capture cross section cr( ca P t ) — to be 
as large as cr( capt ) = 4 x 10" 18 cm 2 . Electron capture 
from residual gas molecules is therefore expected to con- 
tribute significantly to the measured recombination sig- 
nal as a background even at the TSR residual gas pres- 
sure of only 5 x 10~ n mbar. This already proved to be 
all the more the case in a previous TSR experiment with 
Sc 3+ ions Q where due to a technical defect in one of the 
TSR dipole power supplies a reduced maximum rigidity 
of only Bp « 1.2 Tm had been available. Consequently, 
the ion energy had been limited to 300 keV/u. Because 
of the rapid increase of er( capt ) with decreasing ion energy 
and because the experiment had been carried out in early 
summer, i.e., close to a maximum in the slight seasonal 
variation of the TSR residual gas pressure, the measured 
signal to background ratio had been of the order of 1/100. 
Under these conditions the Sc 3+ recombination rate co- 
efficient could only be measured over a limited energy 
range with large statistical uncertainties. In the present 
measurements, that were carried out in December 2000, 
the signal to background ratio was a factor of about 10 
higher as compared to the previous experiment. 



The electron- ion recombination coefficient a(E IC \) is 
obtained as the background-subtracted recombination 
count rate normalized to the electron current and to the 
number of stored ions fl2]| . The estimated systematical 
uncertainties are 15% for the absolute value of the mea- 
sured rate coefficient and less than 2% for the relative 
energy in the energy range under study. 



III. THEORY 

A. Nonperturbative approach of continuum-bound 
transitions 

Starting from the rigorous continuum-bound transition 
theory, Davies and Seaton |2(| discussed the process of 
emission of radiation into the optical continuum due to 
radiative capture of an electron by an ion and derived 
a general formalism including radiation damping. The 
details of their formalism and our numerical scheme for 
application can be found in Ref. |^0| and Ref. [jO], re- 
spectively. Nevertheless, a brief review of the formalism 
and our method follows. 

Using the close-coupling R-matrix approach (see Q 
for details) we start from a treatment of the electron-ion 
collision process that neglects any interaction with the 
radiation field. The resulting wave functions are used 
as a basis for setting up equations including the interac- 
tion with the radiation field in the approximation that 
only the electric-dipole terms are retained and that the 
radiation field is restricted to one-photon and no-photon 
states. The exact solutions for the probability amplitudes 
of the time-dependent matrix equations involved can be 
obtained by the application of a Laplace transform. They 
are expressed in terms of a scattering matrix x with the 
partitioning 



X : 



Xcc Xcp 



(6) 



with \ cc representing the submatrix for electron-electron 
scattering allowing for radiative decays, Xep that for pho- 
toionization, Xpe that for electron capture with the emis- 
sion of a photon, and Xpp that for photon-photon scat- 
tering, and Xpo are written as 



xcc = s[i - 2tt 2 z»(i + zy 1 ^] 

Xpc = -2m{l + Z)- l D\ 



(7) 

(8) 



respectively, where S is the usual electron-electron scat- 
tering matrix in the absence of interaction with radiation 
fields, and D is the reduced dipole matrix (D^ denotes 
its hermitian conjugate) with its matrix element in the 
form 



D 



jJ,j'J' 



2w 3 a 3 
3tt 



1/2 



(lJ II R II l'J') 
(2J+1)V2 



(9) 



where a is the fine-structure constant, lu is the photon 
energy in units of hartrees, and R = rj is the dipole 
operator with the summation extending over all atomic 
electrons. The quantum numbers jJ and 7' J', respec- 
tively, specify the continuum and bound states of the 
atomic system, and J(J') are the total angular momenta. 
The wave function of the continuum electron is normal- 
ized per unit hartree. In Eqs. ^ and || the matrix 



dE- 



(E-tt) 



(10) 
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is related to radiation damping. The variable de- 
notes the total energy of the ion+photon system. In 
the usual first-order theory damping is neglected and 
Xp C = —2ttiD. For the calculation of the PR cross sec- 
tion with damping we emp loy the numerical method de- 
veloped by Zhao et al. [llj] for the evaluation of the prin- 
cipal value of the integral appearing in Eq. [l^. 

In the present calculation the parameters for the Is, 2s, 
2p, 3s, 3p orbitals of scandium were taken from the com- 
pilation of Clementi and Roetti [21j] , and the 3d orbital 
was optimized on the 3p 5 3d l!3 P°, ^ 3 D°, 1 ' 3 F° states 
(weighted) by using the CIV3 code of Hibbert [§|]. We 
also evaluate PR cross sections for further target states 
including 4s, 4p, 4d orbitals. All these orbitals were opti- 
mized in a way similar to the 3d orbital. The symmetries 
2S+i£?r w jtjj jr, < 5 we re included in our calculations. 
Here, S and L are the total spin and orbital angular mo- 
mentum quantum numbers, respectively, and it denotes 
the parity. In our calculation we included configurations 
of odd parity only, since these are expected to lead to the 
dominant resonances in the energy region of our primary 
interest around the 3p 5 3d 2 ( 3 F) 2 F resonance. 



B. Perturbation theory 

On the basis of the principle of detailed balance, the 
photorecombination (PR) cross sections (in atomic units) 
from an initial continuum state j to a final bound state 
/ may be written in the form |23j] 



,T PR 



9f a 2 uJ 2 



9j 



2c 



(11) 



where <?j and <?/ are the statistical weights of the initial 
ion core and the final recombined ion, uj and e are the 
photon energy and the free electron energy in hartrees, 
and cr?? denotes the photoionization (PI) cross section 
from state / to state j. In first order perturbation treat- 
ment it is 



7 PI 



47r 2 aw 



M 



fj\ 



(12) 



where Mfj is the PI matrix element of the corresponding 
transition with the continuum state being normalized per 
unit energy (hartree). According to continuum-bound 
configuration- interaction theory fl24|| , the perturbative PI 
matrix element M_h in low orders may be written in the 
form (see also ||, |2f| ) , 



0W> (!-*£■ 



iTd/2 



E 



(j\V\d)(d\R\f) 
A d + iT d /2 



(13) 



where is the summation of radiative and autoioniza- 
tion widths of the resonance state d, is the autoion- 
ization width from state d to state j, = e — e<i, in 



which ed is the energy level of resonance d, V is the 
electron-electron interaction. It should be noted that in 
the usual independent-processes and isolated-resonance 
approximations all cross terms in |M/., | 2 are omitted that 
lead to interference between resonant and non-resonant 
recombination channels as well as to the interaction be- 
tween resonances. In the present investigation these cross 
terms are taken into account where necessary. 



IV. RESULTS AND DISCUSSION 

A. Experimental results 

Fig. [I] shows the measured Sc 3+ (3s 2 3p 6 *S) recombi- 
nation rate coefficient. The most prominent feature is 
the peaklike structure around E Te \ = 40 eV. It can be 
attributed to mainly unresolved high n Rydberg DR 
resonances attached to 3p 5 3d X P and 3p 5 4s 1,3 P core 
excitations. Vertical lines indicate DR resonance po- 
sitions E{n) as calculated with the Rydberg formula 
E(n) = E[po) - Kq 2 /n 2 with K = 13.606 eV, q = 3, 
and the series limits E(oo) taken from Tab. |. It should 
be noted that we do not observe Rydberg resonances to 
arbitrary high n, since recombined Sc 2+ (3p 6 rii) ions in 
weakly bound high n states are field-ionized by motional 
electric fields in the charge-analyzing dipole magnet and 
therefore do not reach the recombination detector. Ap- 
plying a detailed model of the field-ionizing properties 
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FIG. 1: Measured Sc 3+ (3s 2 3p 6 ) recombination-rate coef- 
ficient. Vertical lines denote resonance positions of the 
3p 5 4s 1 P nl and 3p 5 3d X P nl Rydberg series of DR resonances. 
Vertical arrows denote possible positions of further 3p J 3d se- 
ries limits (see Tab. Q). The insets enlarge the spectrum at low 
energies. Different energy intervals were measured for differ- 
ent time durations and therefore exhibit different statistical 
uncertainties. The rise of the recombination-rate coefficient 
towards E Ic i — eV is due to RR. The full line approaching 
a(E le \) = towards higher energies is the total hydrogenic 
RR rate coefficient (see text). 



5 



of our apparatus that takes into account the radiative 
decay of high Rydberg states on their way from the elec- 
tron cooler to the charge-analyzing magnet as well as 
state-selective field-ionization probabilities Jll| we esti- 
mate the highest detected Rydberg state to be approxi- 
mately n = 20. 

In principle, Rydberg series attached to further 3p 5 3d 
core excitations (Tab. |) could have been expected to 
be visible. The respective core-excitation energies are 
marked by vertical arrows in Fig. |l]. However, no strong 
DR resonances belonging to such series are observed. Ap- 
parently, only dipole allowed (under LS'-coupling con- 
ditions) 1 S^ 1 P core excitations give rise to strong 
DR channels. This has already been observed for iso- 
electronic Ti 4+ ions |^6| where DR proceeds almost exclu- 
sively via resonances attached to 3p 5 3d : P excitations. 
In contrast to the present finding for Sc 3+ , in the DR 
spectrum of Ti 4+ no DR resonances due to 3p— >4s exci- 
tations were found. One strong hint to the origin of this 
difference is the fact that in Sc 3+ the 3p — > 4s excitation 
energies are about 1 eV lower than the 3p 5 3d : P excita- 
tion energy (Tab. ||), while in Ti 4+ they are higher by up 
to 6 eV. Although being iso-electronic Sc 3+ and Ti 4+ do 
have a markedly different structure. 

The recombination rate coefficient rises sharply to- 
wards zero relative energy (lower inset of Fig. Q). This 
rise is due to RR. The RR rate coefficient is rather small 
as can be expected from the fact that RR into the K, L, 
Mi and shells of argonlike Sc 3+ is not possible. This 
is in contrast to e. g. bare Li 3+ ions where these channels 
yield the major contribution (43% at E IC \ — 10~ 6 eV) 
to the total hydrogenic RR cross section |27|] summed up 
to n = 20. The hydrogenic total Sc 3+ rate coefficient, 
i. e. the sum of the nZ-differcntial Li 3+ RR cross sections 
ranging from the 3d-subshcll to the experimental cutoff 



TABLE I: 3p-excitation energies of argonlike Sc 3+ (3s 2 3p 6 ) 
taken from the NIST atomic spectra database J3(| . Note that 
the excitation to the 3p 5 4s terms requires less energy than 
the excitation to the highest 3p 5 3d term. 



Configuration 



Ss^Sp^d 

3s 2 3p 5 3d 

3s 2 3p 5 3d 
3s 2 3p 5 3d 

3s 2 3p 5 3d 
3s 2 3p 5 4s 



3s 2 3p 5 4s 
3s 2 3p 5 3d 



Term 


J 


Energy (eV) 







29.72190 




1 


29.80618 




2 


29.98113 


3 F 


4 


31.08387 




3 


31.25746 




2 


31.41831 


*D 


2 


33.15641 


:i D 


3 


33.23202 




1 


33.40464 




2 


33.40878 


*F 


3 


33.60660 


3p 


2 


41.29801 




1 


41.46096 







41.82486 


!P 


1 


41.84264 




1 


42.77523 



quantum number n — 20 convoluted with the experi- 
mental electron-energy distribution (see e.g. Ref. | 15|) , 
is plotted as the full line in the lower inset of Fig. [lj. 
It agrees with the experimental curve within the experi- 
mental uncertainties. In spite of the long measuring time 
of 60 hours, in particular in the energy range of 8-17 eV, 
the statistical uncertainty did not become low enough to 
bring out all structures. At the present level of statis- 
tical uncertainty it cannot be decided whether the faint 
structures in the energy range 1-7.5 eV are due to weak 
DR resonances or due to unaccounted background mod- 
ulations. 

The DR resonances visible in the upper inset of Fig. |l| 
are due to 3p 5 3d 2 and 3p 5 3d4/ doubly excited configura- 
tions. These configurations are considerably shifted down 
in energy as compared to the Rydberg energies shown by 
vertical lines in Fig. [l] and the energies of the various 
terms must be taken from detailed calculations. As al- 
ready mentioned in the introduction the terms belonging 
to the 3p 5 3d 2 configuration straddle an energy range of 
about 20 eV. Also the 3p 5 3d 4/ configurations can be ex- 
pected to be subject to strong fine structure interactions 
and hence to be spread out over several electron volts. 
The energy range between about 8 and 26 eV (upper in- 
set in Fig. [l]) was searched preferentially in order to de- 
tect the broad 3p 5 3d 2 2 F Super-Coster-Kronig resonance 
predicted by theory. The measurement reproduces the 
strong double peak structure near 15 eV found already 
in the previous experiment JtJ (scan range 12-18 eV) but 
now also reveals considerable additional structure both 
above and below these peaks. While much of this struc- 
ture seems to arise from a superposition of several nar- 
row peaks, a single broad resonance (with some narrower 
structure superimposed on its low energy slope) seems to 
dominate the recombination rate near 12.3 eV. This sig- 
nal bears strong similarity with the predicted 3p 5 3d 2 2 F 
Super-Coster-Kronig resonance || and we take this fea- 
ture as a key to assigning a considerable part of the ob- 
served spectrum. This interpretation is also supported by 
the argument that for the isoelectronic system Ti 4+ |2(| 
the 3p 5 3d 2 2 F resonance was observed about 3 eV below 
the prediction obtained by the same theoretical approach 
|p8| , so here a position of ~ 12 eV can be expected. 

A closeup of the experimental spectrum is shown in 
Fig. |^, where the resonances in the energy range 10.9- 
16.5 eV have been fitted by either Lorentzian line profiles 



A 2 

L(E) = —-- 

T d 7T 1 



or an asymmetric Fano profile 



F(E) = 



A 2 



Q 2 r d 7r 



(Q + sf 



(14) 



(15) 



with e = 2(E—Ed) /Td- The term -1 in the square brackets 
ensures that F(E) — > for E — ► ±oo and the normaliza- 
tion factors have been chosen such that F(E) — > L(E) 
in the limit Q — > oo. Conversely, small absolute val- 
ues of the asymmetry parameter Q lead to significantly 
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FIG. 2: Fit of Lorentzian and Fano line profiles to the exper- 
imental data over the energy range 10. 9-16. 5 eV. The individ- 
ual line profiles have been convoluted with a Gaussian repre- 
senting the experimental energy distribution (see appendix). 
The fit results are summarized in Tab. [n|. The fit delivers a 
significant asymmetry (Q = 6.3±1.8)of the broad resonance 
at 12.3 eV. Moreover, from the fit a longitudinal electron tem- 
perature ksTu = 0.15 ± 0.03 meV is obtained. 



asymmetric line shapes. Assuming that the dominant 
contribution to the resonance width stems from au- 
toionization and neglecting contributions from neighbor- 
ing resonances, Q can be expressed in terms of the matrix 
elements occurring in Eq. fL3l i. c., 



Q = 



2(j\V\d) (d\R\f) 
r d (j\R\f) 



(16) 



as can be shown by simple algebraic manipulations on 
|M/j| 2 . For each resonance the parameters varied in the 
fit were the peak area A, the resonance energy Ed, the 
resonance width Yd, and in case of the 3p 5 3d 2 2 F res- 
onance the asymmetry parameter Q. In the fit both 



TABLE II: Results obtained from fitting Lorentzian (Eq. |l^) 
and Fano (Eq. [IB]) line profiles (convoluted with a Gaussian, 
see appendix) to resonances in the energy range 10.9-16.5 eV 
(see Fig. ^). The errors given do not include systematic un- 
certainties in the energy calibration (±2%) and the absolute 
rate coefficient determination (±15%). 



E d 
(eV) 




r 

(eV) 




(10 


-12 


A 

cm 3 s _1 eV) 


11.175 ± 


0.041 


0.317 ± 


0.099 


3.8 


± 


1.8 


11.481 ± 


0.143 


0.140 ± 


0.253 


2.6 


± 


6.1 


11.587 ± 


0.036 


0.000 ± 


0.193 


1.2 


± 


3.8 


11.819 ± 


0.004 


0.000 ± 


0.021 


1.9 


± 


0.4 


12.305 ± 


0.027 


0.889 ± 


0.066 


14.3 


± 


1.0 a 


14.802 ± 


0.001 


0.056 ± 


0.006 


11.6 


± 


0.2 


15.324 ± 


0.001 


0.018 ± 


0.006 


9.1 


± 


0.2 



"Fano profile (Eq. |l|) with Q = 6.3 ± 1.8. 



L(E) and F{E) have been convoluted (see appendix) 
by a Gaussian representing the experimental energy- 
distribution function with its full width at half maximum 
(FWHM) AE re A(E d k-BT\\ hx2) 1 ' 2 in the energy region 
of interest jl^] . The longitudinal electron-beam tempera- 
ture /cbTh (&B denotes the Boltzmann constant) has also 
been allowed to vary during the fit. From the fit we 
obtain k-gT^ = 0.16 ± 0.03 meV. Further fit results are 
summarized in Tab. |n|. 

The most important result of the fit is the determi- 
nation of the 3p 5 3d 2 2 F resonance energy and width 
to 12.31 ± 0.03 eV (±2% systematic uncertainty) and 
0.89 ± 0.07 eV, respectively. For the asymmetry param- 
eter of this resonance we obtain Q = 6.3 ± 1.8. A ma- 
jor source of uncertainty beyond the errors given is the 
background from collisions with residual gas molecules. 
Under the conditions of the present experiment the back- 
ground causes the major fraction of the measured recom- 
bination count rate (see Sec. |l]). A subtraction of this 
high background level from the measured signal is dif- 
ficult, especially, if changes of the residual gas pressure 
have occurred on the time scale of the switching between 
signal and background measurement. In this case a small 
fraction of the background may not have been properly 
subtracted possibly leading, e. g., to the nonzero recombi- 
nation signal in the energy range 15.6-17 eV. In principle, 
this signal might also be caused by a group of weak DR- 
resonances, however, this assumption is not supported by 
our theoretical calculations (see below) . When allowing a 
constant background level as an additional free parame- 
ter in the fit a considerably different asymmetry parame- 
ter Q = 10.4 ±4.1 is obtained with the fitted background 
level amounting to 6.0 ± 0.9 x 10 -13 cm 3 s _1 . On the 
other hand the resonance position and width change only 
within their errors as listed in Tab. |TJ to 12.35 ± 0.03 eV 
and 0.85 ± 0.07, respectively. 

The significance of our fit result for the asymmetry pa- 
rameter may also be assessed by fitting the 3p 5 3d 2 2 F res- 
onance with a symmetric Lorentzian instead of an asym- 
metric Fano line shape. From such a Lorentzian fit we 
obtain a x 2 value that is higher by only 5% as compared 
to the Fano fit, i.e., the experimental evidence for an 
asymmetric line shape is not very strong. 

Fig. U shows the comparison of the present experimen- 
tal result with the theoretical result of Gorczyca et al. 
H in the region of the 3p 5 3d 2 2 F resonance. As already 
inferred from our previous work the theoretical reso- 
nance positions do not match the experimental ones. For 
the comparison presented in Fig. |3| the theoretical en- 
ergy scale has been shifted by 2.2 eV in order to match 
the experimental position of the broad 3p 5 3d 2 2 F res- 
onance. Moreover, the absolute theoretical rate coeffi- 
cient has been normalized to the experimental height of 
this resonance by a multiplication factor of 0.63. After 
these adjustments a basic similarity between predicted 
and measured spectrum is found. In particular, a broad 
feature is observed strongly resembling the calculated 
3p 5 3d 2 2 F resonance. However, in the experimental spec- 
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FIG. 3: Comparison of the measured Sc 3+ rate coefficient 
with the theoretical result of Gorczyca et al. |g] (full curve) in 
the region of the 3p 5 3d 2 2 F resonance. The theoretical peak 
designations are given. The theoretical rate coefficient has 
been shifted by 2.2 eV to lower energies and multiplied by 0.63 
in order to obtain agreement with the measured 3p 5 3d 2 2 F 
resonance position and height. 



trum considerably more resonances are visible than in the 
theoretical one. The theoretically predicted single peaks 
at 11.6 and 15.1 eV appear to be split into two compo- 
nents. Also, the predicted relative resonance strengths 
deviate from the experimental findings. The theoretical 
3p 5 3d 2 2 P peak at 15.1 eV is much smaller than the ex- 
perimental doublet. 



B. Theoretical results 

The result of our LS-coupling-R-Matrix calculation of 
the total Sc 3+ PR cross section including damping is plot- 
ted in Fig. [|a. It can be noticed that there is a very wide 
resonance at about 13.6 eV. As already predicted by Gor- 



TABLE III: Comparison of experimental and theoretical en- 
ergy levels of some doubly excited Sc 2+ states relative to the 
3p 6 1 S state in units of eV. The assignment of the experimen- 
tal resonances below 12 eV is uncertain. 



State 


Exp. 




Theory" 








A 


B 


C 


D 


3p b 3d( 3 F)4s 2 F " 

3p 5 3d( 1 D)4s 2 D 

3p 5 3d( 1 F)4s 2 F 

3p 5 3d( 3 D)4s 2 D 

3p 5 3d 2 ( 3 F) 2 F 

3p 5 3d 2 ( 3 P) 2 P 

3p 5 3d 2 ( 3 F) 2 D 


11.18 
11.48 
11.58 
11.82 
12.31 


11.00 
12.71 
13.61 


14.62 


11.59 
12.18 
12.77 
12.89 
13.23 


15.74 


14.80 
15.32 


16.40 


17.35 


15.85 
16.90 


16.81 




10 15 20 

Electron energy (eV) 




5 10 15 20 25 

Relative energy (eV) 



FIG. 4: (a) Total Sc 3+ photorecombination cross section 
as calculated with the R-matrix approach in the LS cou- 
pling scheme including 3p 5 3dnZ DR channels with n < 4 
and with the symmetries 2 L n with L < 5. (b) Merged 
beams recombination-rate coefficient derived from the above 
R-Matrix cross section by convolution with the experimental 
electron-energy distribution function with ksT^ = 0.15 meV 
and k B T ± = 10 meV. 



czyca et al. @, it is due to DR via the 3p 5 3d 2 ( 3 F) 2 F 
state. However, our resonance energy is by about 1 eV 
closer to the experimental value of 12.29 eV (see Tab. ||) 
than that of Gorczyca et al. This difference originates 
from the use of different sets of basis functions in both 
calculations. 



"A: this work, B: Ref. ||, C: Ref. S, D: Ref. 
i, ls 2 2s 2 2p 6 3s 2 is omitted. 



Tabic III lists measured and theoretical 
resonance energies including further theoretical results 
pTl|. [Io[ . The relatively large differences between the var- 
ious theoretical results illustrates again the difficulty to 
exactly describe the many-electron atomic system under 
consideration. 

In order to be able to more accurately compare the 
theoretical result with the measured merged-beams rate 
coefficient we convoluted the R-matrix cross section with 
the experimental electron energy distribution function, 
that can be parameterized by the longitudinal and trans- 
verse (with respect to the electron beam direction) elec- 
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FIG. 5: Comparison of the measured Sc 3+ rate coefficient 
with the present theoretical result (R-Matrix + perturbation 
theory) in the region of the 3p 5 3d 2 ( 3 F) 2 F resonance. The 
theoretical peak designations are given. The theoretical rate 
coefficient has been shifted by 1 .13 eV to lower energies and 
multiplied by 0.55 in order to obtain agreement with the mea- 
sured 3p 5 3d 2 ( 3 F) 2 F resonance position and height. 



tron beam temperatures Tu and T±, respectively [ [T2| . 
For the convolution result displayed in Fig. |]b we have 
used fc B 7|| = 0.15 meV (cf., Fig. |) and fc B T ± = 10 meV. 
Comparing with the experimental data in Fig. [l] we find 
that the R-matrix calculation does not reproduce the ex- 
periment in every detail. For example, the calculated 
resonance peaks in the energy range 7-17 eV are too 
high, whereas those beyond 19 eV are too low. More- 
over, some of the observed peaks at lower energies are not 
reproduced by the R-matrix calculation. The disagree- 
ment at high energies can be traced back to the neglect 
of even-parity configurations and to the neglect of rela- 
tivistic effects by restricting the R-Matrix calculation to 
LiS-coupling. This also explains the missing of peaks at 
lower energies. Using Cowan's code || for the calcula- 
tion of DR r esona nces in the framework of perturbation 
theory (Sec. Ill B| ) we find that due to relativistic effects 
(mainly spin-orbit interaction) also LS'-forbidden excita- 
tions to 3p 5 3d 2 ( 3 F) 2 D as well as to 3p 5 3d( 1 D)4s 2 D and 
3p 5 3d( 3 D)4s 2 D doubly excited states may contribute to 
the PR cross section. It should be emphasized that res- 
onances arising from relativistic effects have in general 
small Auger widths if the atomic number Z is not large. 
As has been explained in some detail by Zhao and Shirai 
pl[ they, nevertheless, may still play an important role 
in recombination if their radiative lifetime is of the order 
of the radiative lifetimes of the dominant DR resonances. 

As shown in Fig. |S[ a much better agreement with 
the experimental data is achieved when the the 2 D reso- 
nances mediated by spin-orbit effects are taken into ac- 
count in addition to the R-Matrix cross section. For 
the generation of the incoherent sum of the R-Matrix 
+ perturbation-theory rate coefficient shown in Fig. |5| 




11 13 15 

Relative energy (eV) 



FIG. 6: (a) Cross section for recombination of Sc via the 
intermediate 3p 5 3d 2 ( 3 F) 2 F state into the 3p 6 3d 2 D state: to- 
tal cross section (solid line), nonresonant RR cross section 
(dashed line), DR cross sections (dotted line), and interfer- 
ence term (dashed-dotted line) . Also given are the total cross 
section and the interference term for the case of an artificially 
reduced (by a factor 10) RR cross section (thinner lines), (b) 
Merged-beams recombination rate coefficients corresponding 
to the calculated (dotted line) and artificial (solid line) total 
cross sections. 



we have used the 2 D Auger and radiative rates as cal- 
culated by Cowan's code. The resonance energy, how- 
ever, has been determined by using the experimentally 
observed energy splittings between the 2 D resonances and 
neighboring resonances obtained from R-Matrix theory. 
For example, the experiment yields 3p 5 3d 2 ( 3 P) 2 P and 
3p 5 3d 2 ( 3 F) 2 D resonance energies of 14.8 and 15.3 eV, 
respectively (Tab. [n]). The level difference equals 0.5 eV. 
Since our R-matrix calculation yields a 3p 5 3d 2 ( 3 P) 2 P 
resonance energy of 16.4 eV, we chose 16.9 eV to be the 
3p 5 3d 2 ( 3 F) 2 D resonance energy. Moreover, for the com- 
parison in Fig. H the theoretical merged-beams rate co- 
efficient was multiplied by a factor 0.55 and shifted by 
1.13 eV to lower energies. The difference between exper- 
imental and theoretical cross sections of nearly a factor 2 
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FIG. 7: Partial cross sections for recombination via 2 F° 
states [in the order of increasing energy: 3p 5 3d( 3 F)4s 2 F, 
3p 5 3d( 1 F)4s 2 F and 3p 5 3d 2 ( 3 F) 2 F] calculated in the R-matrix 
approach, showing the variation of the resonance line shapes 
as a function of energy separation (a) as calculated ab initio, 
(b) 0.5 times, and (c) 1.5 times the separation in (a) with no 
other parameter varied. 



is found to be independent of the inclusion of 4s, 4p and 
4d orbitals in the set of target states. 



C. Line shape of the 3p 5 3d 2 ( 3 F) 2 F resonance 

Clearly, the theoretical line shape of the broad 
3p 5 3d 2 ( 3 F) 2 F resonance is asymmetric. As already men- 
tioned in the introduction, Gorczyca et al. H have 
attributed this asymmetry to interference between RR 
and DR pathways. However, this interpretation may be 
doubted in view of the fact that in the relevant energy 
range the RR cross section is orders of magnitude smaller 
than the DR cross section. Naively, one would expect 
that a significant interference between DR and RR re- 
quires cross sections of similar magnitude. It must be 
kept in mind, however, that the asymmetry parameter 



Q is also governed by the natural line width Td- Inspec- 
tion of Eq. [l6] shows that a large value of T<i can make 
up for a small RR transition matrix element (j\R\f). 
Fig. ||a shows the RR, DR and interference contributions 
to the recombination cross section in the energy region of 
the 3p 5 3d 2 ( 3 F) 2 F resonance as calculated by perturba- 
tion theory. With the width of the 3p 5 3d 2 ( 3 F) 2 F reso- 
nance due to its Super-Coster-Kronig decay channel be- 
ing extraordinarily large (Fa = 1.06 eV obtained from 
Cowan's atomic-structure code [|| ) the asymmetry in the 
merged-beams rate coefficient is retained even when the 
RR cross section is artificially lowered by a factor of 10 
(Fig. ||b). The theoretical asymmetry parameter changes 
from Q — 3.97 to Q — 12.6. Our experimental findings 
of T d = 0.89 eV and of Q = 6.3 (see Tab. ||) correspond 
to these theoretical results. 

Apart from interference between DR and RR another 
source of distorted line shapes is the interaction be- 
tween overlapping DR resonances of the same symmetry 
pS2f . From our R-Matrix calculation we find two fur- 
ther resonances of 2 F symmetry, namely the resonances 
3p 5 3d( 3 F)4s 2 F and 3p 5 3d( 1 F)4s 2 F located at 11.00 and 
12.71 eV, respectively. As shown in Fig. Ra, especially the 
latter is overlapping with the broad 3p^3d 2 ( 3 F) 2 F reso- 
nance. In Figs. [7b and |?]c we show how the line shape 
of the 3p 5 3d 2 ( 3 F) 2 F resonance depends on its (artifi- 
cially varied) distance to the 3p 5 3d( 1 F)4s 2 F resonance. 
Clearly the shapes of both resonances are affected by 
their mutual separation. The experimental observation 
of this effect, however, is hindered by the presence of 
the 2 D resonances close to the 3p 5 3d( 1 F)4s 2 F resonance. 
Under this circumstance the experimental resolution does 
presently not suffice to reveal the theoretically predicted 
subtle resonance-resonance interference. 



V. CONCLUSIONS 

In the experimental energy range 0-45 eV photore- 
combination of Sc 3+ (3s 2 3p 6 ) ions is dominated by the 
3p 5 3d 1 Pn^ and 3p 5 4s 1 Pn/ Rydberg series of DR res- 
onances. Significant differences to the recombina- 
tion of isoelectronic Ti 4+ ions (26| are found that are 
due to subtle differences in the electronic structure of 
Sc 3+ and Ti 4+ . Level energies, widths and resonance 
strengths were determined for a number of doubly excited 
Sc 2+ (3s 2 3p 5 3d 2 ) and Sc 2+ (3s 2 3p 5 3d4s) states. Further- 
more, we find a small experimental evidence for an asym- 
metry of the 3p 5 3d 2 ( 3 F) 2 F resonance line shape. A 
more precise determination of the line shape from the 
experimental data would require less statistical uncer- 
tainties especially in the high energy tail of the reso- 
nance. In the present measurement the high background 
level from electron capture in Sc 3+ collisions with resid- 
ual gas molecules prevented this requirement from being 
fulfilled. Moreover, the peaks on the low energy side of 
the 3p 5 3d 2 2 F resonance leave additional room for ambi- 
guity. 
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Starting from rigorous continuum-bound transition 
theory the Sc 3+ photorecombination cross section has 
been calculated using both the R-matrix method and a 
perturbative treatment. We found that the line shape 
of the 3p 5 3d 2 ( 3 F) 2 F resonance is influenced by both in- 
teracting resonances and interference between RR and 
DR. In spite of the fact that the Sc 3+ RR cross section is 
almost negligibly small, DR-RR interference becomes no- 
ticeable due to the large 3p 5 3d 2 ( 3 F) 2 F resonance width 
of 0.89 eV. Moreover, we have shown that relativistic 
effects play an important role in the recombination of 
Sc 3+ . Three experimentally observed strong 3p 5 3d 2 2 D 
and 3p 5 3d4s 2 D resonances can only be reproduced the- 
oretically by invoking relativistic effects. Our R-matrix 
calculations as well as the calculations of Gorczyca et at 
H that both have been restricted to LS-coupling fail in 
reproducing these strong experimental peaks. 

Our results show the limits of the independent- 
processes and isolated-resonances approximations usu- 
ally made in photorecombination calculations. A more 
stringent test of our theoretical results could in principle 
be provided by the measurement of high-resolution pho- 
toionization cross sections of Sc 2+ (3p 6 3d 2 D) ions. Re- 
sults of such an experiment will be reported elsewhere 
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APPENDIX A: CONVOLUTION OF A 
FANO-PROFILE WITH A GAUSSIAN 

The convolution of a Fano profile (Eq. |l5|) with a nor- 
malized Gaussian 



4(ln2)£ 2 
(AE) 2 



(Al) 



with AE being its full width at half maximum (FWHM), 
yields 



/oo 
F(E') G{E' - E) dE' 
-OO 

[Qy + (t- x)] 2 

, (t-x) 2 + y 2 




(A2) 
e'^dt- 1 



AE^H 
where the definitions 



1-^W(*)]-|9K*)] 



2 vines' - E) 
AE ' 
2V\n2(E d - E) 



V 

z 




(A3) 

(A4) 

(A5) 
(A6) 



have been used and where w(z) denotes the complex error 
function (34) . Its real and imaginary parts can efficiently 
be calculated with a numerical algorithm by Humlicek 
fl35[ . It should be noted that for Q — > oo Eq. ^2| yields the 
Voigt profile, i. e. the convolution of a Lorentzian (Eq. |l4| ) 
with a Gaussian. 
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